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Abstract 

A non-universal scalar mass supergravity type of model is explored where the first 
two generation of scalars and the third generation of sleptons may be very massive. 
The lighter or vanishing third generation of squarks as well as Higgs scalars at the 
unification scale cause the radiative electroweak symmetry breaking constraint to be 
less prohibitive. Thus, both FCNC/CP-violation problems as well as the naturalness 
problem are within control. We identify a large slepton mass effect in the RGE of ni^^ 
(for the down type of Higgs) that may turn the later negative at the electroweak scale 
even for a small tan j3. A hyperbolic branch/focus point like effect is found for m\ that 
may result in very light Higgs spectra. The lightest stable particle is dominantly a bino 
that pair annihilates via Higgs exchange, giving rise to a WMAP satisfied relic density 
region for all tan (3. Detection prospects of such LSPs in the upcoming dark matter 
experiments both of direct and indirect types (photon flux) are interesting. The Higgs 
bosons and the third generation of squarks are light in this scenario and these may be 
easily probed besides charginos and neutralinos in the early runs of LHC. 
PACS No: 04.65. +e, 13.40Em, 14.60Ef, 13.85.-t, 14.80.Ly 

1 Introduction 

Low energy supersymmetry (SUSY) [1] is one of the most active fields of research for physics 
beyond the standard model (SM) [2]. A minimal extension of the Standard Model when 
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super symmetry is incorporated is the Minimal Supersymmetric Standard Model (MSSM) 
[1,3] that includes two Higgs doublets. The model however has a large number of SUSY 
breaking parameters and this motivates one into studying models with specific mechanisms 
for breaking SUSY. The later involves high scale physics input and renormalization group 
analyses. This in general leads to a large reduction of the number of unknown parameters. 
The minimal supergravity (mSUGRA) [4] model is one of the well studied SUSY models. It 
requires a very few input parameters at the gauge coupling unification scale or grand unified 
theory (GUT) scale, Mq — 2x10^^ GeV. The parameters indeed quantify our ignorance of the 
exact nature of SUSY breaking. The model incorporates radiative breaking of electroweak 
symmetry (REWSB). The unification scale universal input parameters are: i) the gaugino 
mass parameter mi, ii) the scalar mass parameter mo and iii) the tri-linear SUSY breaking 
parameter Aq. Additionally, one has to provide with tan/5, the ratio of Higgs vacuum 
expectation values and the sign of the Higgsino mixing parameter fi. Renormalization group 
evolutions are used to obtain the electroweak scale parameters of MSSM. 

There is however no a priori necessity of considering such universalities of parameters. 
Indeed, it is worthwhile to explore scenarios with non-universalities in the scalar or in the 
gaugino masses at the unification scale [5-11] . Non-universal scalar masses may appear 
because of non-fiat Kahler potential [12]. However, one must be careful to accommodate the 
stringent constraints from phenomena involving fiavor changing neutral currents (FCNC). 
Satisfying FCNC constraints demands near-degeneracy of the first two generations of scalar 
masses but the requirements on the third generation of scalars as well as the Higgs scalars are 
not so stringent [13,14]. Thus it is seen that FCNC constraints as well as constraints from CP- 
violating phases (for example, those arising from the electric dipole moments of the electron 
and neutron) may be managed by introducing multi-TeV scalar masses for the first two 
generations of scalars [15]. The third generations of scalar masses and the Higgs scalar masses 
however should be adequately light in order to satisfy naturalness [16]. There have been 
several efforts for obtaining the desired features as outlined above. Analyses with Radiatively 
Generated Inverted Mass Hierarchy Models (RIMH) were made in Ref. [14,17,18]. In Ref. [18] 
the authors achieved the above-mentioned requirements at the electroweak scale by using 
t—b—T Yukawa unification with special non-universal relationship among the scalar masses at 
Mq. Here, the Higgs and the third generation of scalar masses are rapidly diminished at the 
electroweak scale via RG evolutions. However, the Yukawa unification and the consideration 
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of REWSB constrain such models heavily. A second realization of the above idea is partially 
possible via the hyperbolic branch (HB)/ focus point (FP) [19,20] scenarios where the scalars 
may become considerably massive (multi-TeV) in a subset of the typical mSUGRA parameter 
space satisfying universal boundary conditions while fine-tuning [16] still remains small. A 
third possibility was considered in Ref. [21], where plain decoupling arguments motivated 
the authors in using explicit splitting at Mq between the scalars belonging to the first two 
generations and the same of the third generation (along with the Higgs scalars). Here the first 
two generations of scalars were chosen in the multi-TeV domain whereas the third generation 
of scalars as well as the Higgs scalars were considered to be in the sub-TeV zone. Additionally, 
a large value of the tri-linear coupling parameter was chosen in Ref. [21] for the first two- 
generations. Universality of scalar masses in the first two-generations along with a choice 
of a different scalar mass parameter for the third generation as well as the Higgs scalars, 
or even splitting of squarks and sleptons within the third generation itself have also been 
considered in Refs. [22-24] . In Ref. [23] the authors additionally considered non-universality 
in the gaugino masses and analyzed the fine-tuning aspect of computing the relic density 
of dark matter in addition to obtaining a parameter zone of the MSSM corresponding to a 
well tempered neutralino [25]. Similarly, analyses with only non- universalities in the Higgs 
scalar may be seen in Refs. [6,8,10,26]. A comprehensive set of characteristic possibilities 
for varieties of non- universal SUGRA scenarios may be seen in Ref. [27]. 

In this analysis we explore a SUGRA scenario with non-universal scalar masses that i) 
allows to have very large first two-generation of scalar masses so as to obey the FCNC and 
the CP-violation limits easily {ie. without requiring any ultra-small phases) and that would 
not impose any additional price on fine-tuning, ii) spans a large amount of MSSM parameter 
space satisfying the neutralino relic density constraint from WMAP data by not requiring 
any delicate mixing of bino and Higgsinos, so that we would find a bino-dominated lightest 
neutralino for most of the parameter space and, iii) satisfies the Higgs mass lower bound from 
LEP2 data as well as other low energy constraints. Certainly, a model with REWSB and 
universal scalar mass like mSUGRA is not friendly to achieve these objectives if we consider 
a common gaugino mass parameter below a TeV or so. The above phenomenologically 
inspired requirements in combination motivate us to introduce non-universality between the 
third and the first two-generation of scalars as well as the Higgs scalars. This has to be 
such that the REWSB conditions would not become prohibitive to have a multi-TeV first 
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two-generation of scalars. Here we would like to point out that we would not be able to 
satisfy the mentioned objectives by considering non- universalities only in the Higgs scalar 
masses. 

A simpler and purely phenomenological attempt in this direction in a universal gaugino 
mass framework could be to generate all the third generation of scalar masses and the Higgs 
scalar masses radiatively, starting from zero values at the unification scale, while keeping the 
first two-generations of scalar masses universal (and this may be in the multi-TeV zone). We 
point out that starting with a similar range of values for the Higgs scalars as with the third 
generation of scalars at Mq would be desirable since this would lead to similar range of mass 
values for all the soft-SUSY breaking terms contributing to REWSB after RG evolutions. 
However, if we consider all the third generation of scalars at Mq to be light, we would find 
stau (f) becoming the lightest stable particle (LSP) or even tachyonic at the electroweak 
scale. On a similar note, we remind that such an appearance of tachyonic sleptons also arise 
in the Anomaly Mediated Supersymmetry Breaking (AMSB) analyses, and it is avoided 
purely by a phenomenologically inspired way of adding an appropriate non- zero mass value 
to all the scalars at the unification scale in the minimal AMSB scenario [28]. Thus with a 
simple motivation of managing with FCNC and CP-violation, naturalness as well as the dark 
matter constraints concurrently we consider a non-zero mass value for the third generation 
of sleptons, while having the masses of squarks of the third generation and the Higgs scalars 
vanishing at Mq- Additionally, for convenience we set the third generation of slepton mass 
parameters the same as that of the first two-generation of scalars at Mq- 

Thus, the parameters of our Non-universal Scalar Mass model which will henceforth be 
called the NUSM model is given by 



where the scalar mass input assignments at Mq are as follows. 

i) The unification scale mass parameter for the first two-generations of squarks, sleptons and 
the third generation of sleptons is mo, where mg is allowed to span up to a very large value. 

ii) The mass parameters for the third generation of squarks and Higgs scalars are set to zero. 
We could have also chosen a non- vanishing value for ii) as long as it is sufficiently small. We 
may note here that different values of mass parameters for squarks and sleptons at Mq niay 
appear in orbifold models with large threshold corrections. This is related to having different 
modular weights associated with the squarks and the sleptons in a given generation [29,30]. 



tan/3, mi , Aq, S2(?n(/i), and mo 
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As we will see below, the NUSM model provides with a highly bino-dominated neutralino 
dark matter over almost its full range of parameter space. It has an interesting feature of 
having a large funnel region, a region of parameter space where the associated annihilation 
channels are characterized by the direct channel pole 2171^0 ~ m^, tuh- We note that unlike 
mSUGRA where one finds the funnel region only for a large value of tan f], here in the NUSM 
one finds it for almost all possible values of tan/?(^ 5). We will also see that the NUSM 
is further characterized by a lighter or niH particularly when ttiq is large and this is 
found even for small values of tan/3(~ 10). A small part of the parameter space may be far 
from the decoupling [31] region of Higgs boson mixing that causes the lighter CP-even Higgs 
boson (/i-boson) to be non-Standard Model like [32]. This in turn reduces the lower bound 
of rrih much below the LEP2 Higgs boson mass limit. 

In this work we include a semi-analytic calculation that first points out a large mass 
effect in the solution of the renormalization group equations (RGE) of fn^Ho- "^^^ effect 
causes a hyperbolic branch/focus point like behavior^ in m^(~ tti^Ho ^^at may cause 

ttla to become smaller even for a small tan /3 as mentioned above. Here we remind that ttla 
typically becomes smaller in mSUGRA only for large tan (3. 

The paper is organized as follows. In Sec. 2 we will primarily discuss the large mass RGE 
effects in the NUSM on m\ and its consequent reduction for large mo via an HB/FP-like 
effect. We will obtain the semi-analytic results which will also be verified by numerical 
computation. In Sec. 3 we will study the cold dark matter constraint from WMAP data 
including also the constraints from the LEP2 Higgs bound, 6 — > s + 7, and Bs — > We 
will also discuss a few sample points in the context of LHC reach. In Sec. 4 we will discuss 
the direct and indirect detection (via continuous gamma ray) rates of the LSP. Finally we 
will conclude in Sec. 5. 

2 The non-universal scalar scenario: NUSM 

In this section we focus on the key elements that are important for the NUSM. The primary 
quantities of phenomenological interest are /i, and tua- These in turn have important signif- 

^ We remind that the wcU known HB /FP effect that occurs in m^j^ is associated with the first minimization 
condition of Radiative electroweak symmetry breaking. In this analysis we do not have such an effect. On the 
contrary we have a similar HB /FP effect for small values of tan (3 in connection with the second minimization 
condition that is associated with m\. 
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icance on the following: i) the issue of fine-tuning, ii) the presence of non-Standard Model 
like lighter Higgs boson mass bound for a limited region of parameter space, iii) dark matter 
and iv) collider discovery possibilities. Additionally, we remind ourselves about obtaining a 
lighter third generation of squarks at the electroweak scale and this is of course related to 
their vanishing values at the unification scale. 

To start with we write down the REWSB results, 

2_ 1,^2 , rnl^-ml^tan^ /3 ^ Si-Sstan^/? 
^ --^^^^ tan2/5-l + tan' (3-1 ' 



and, 

sin 2p = 2Bfx/{ml^ + m^^ + 2fx' + Si + S2) (3) 

where Si 2 represent the one-loop corrections [33,34], that become small in the scale where 
the Higgs potential Vniggs is minimized. We may approximately consider /i^ ^ —m^j^ (for 
tan /5 ^ 5) and m\ = Tn\j^ + m^u^ + ^ '^'hd "~ ^^Hu ^^^^ level. 

The associated one-loop RGEs are given by (neglecting the very small first two-generations 
of Yukawa contributions), 

dm? 3 

— = (3a2m^ + -aifhl) - 3Yh{mjj^ + m^^ + m% + A^'^) 

3 

-Y^{mj^^ + ml + ml + A^') + —aiSo (4) 



^^^^^ = (3a2m^ + ^aiml) - 3Yt{mjj^ + m^ + m^ + At') - ^( 



Here we have t = ln{MQ/Q') with Q being the renormalization scale. = aj/(47r) for 



z = 1, 2, 3 are the scaled gauge coupling constants (with ai = |ay) and rhi are the running 
gaugino masses, where z = 1, 2, 3 refers to U(l), SU(2) and SU(3) gauge groups respectively. 
Yj represent the scaled and squared Yukawa couplings, e.g, Yj = h'/[4:T\:)' where hj is a 
Yukawa coupling (j = 1, 2, 3 stands for t, 6, r respectively). The quantity Sq which is shown 
below may become relevant for a general set of non-universal boundary condition for scalars. 
However in this analysis it is zero at Mq and it causes only a negligible effect. 

'S'o = - mjj^ + Sfc (mj^ + m^^ + mj^ - m^ - 2m\ ) + (m^ + + m| - - 2m^) (6) 

where the subscript k for /c = 1 or 2 indicates the first two generations. 
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The solutions for /i^ and valid up to a moderately large value of tan/3 (10 or so) follow 
from Ref. [6] after appropriately considering the NUSM parameters. This involves ignoring 
hh and h-r in the limit of a small tan/3. The solutions read 

IJ,^ = mlCi + AlC2 + m\C3 + miAoC4-^Ml + ^^^J^^Sop, and (7) 

1 6 tan^ 3 + 1 

m\ = mlDi + AID2 + mlD^ + miA^D^ - -Ml + -—^ — -Sop. (8) 

22 Z o tan p — 1 

Here p is given by: p = ^[1 — (^^(Sy)]- '^^^ unification scale conditions for scalar mass 
squares for the NUSM may be given as mf = (1 + 5i)irn^. Here, the subscript i for 5i stands 
for: i = Hd, Hu, Ql, ur, dn, II, cr. For NUSM one has 6h^ = = ^Ql = = = -1 
and 61 = (5ejj = 0. These unification scale conditions allow us to define a quantity 6 that 
appears in the expressions of the coefficients Ci and Di as shown below. It follows that 
S = -1 for the NUSM, and 5 = for mSUGRA. We note that only Ci among d's and Di 
among Dj's depend on 6 and one obtains^. 



1 on 1 

Ci = (1 + ^) , 2n J l-^M^tan^/3), and 
tan^ p — 1 2 



Here, — 1 — (m(/200 sin/?)^ ^ 0.2. The definitions of the quantities Ci and Di for i 1 
may be seen in the appendix. As mentioned above, for the NUSM one has 6 = —1. Hence, 
Eqs. 7-9 show that at the level of approximation where hb and hr are ignored, /i^ and m\ 
are independent of rriQ. We will however see that both /i^ and m\ would actually depend on 
rriQ. While the dependence of //^ on tuq would be caused by two-loop RGE effects in rn'jj^, 
the same for m\ is very prominent even at the one-loop level of RGE and this is found when 
we include the effects due to hi, and hr which are specially important for the NUSM. The 
following subsection describes our improved results. 



■^In Ref. [6] the authors considered non-universalities in the masses of Higgs scalars and the third gener- 
ation scalars mqj^ and m„jj. However, NUSM additionally requires non- universality in m^j,. 
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2.1 Large slepton mass RGE effect on pseudoscalar Higgs boson— 
A hyperbolic branch/focus point like effect in m\ for small 

tan P 

Numerical computation shows that unhke what is seen from Eqs. 8 and 9, m\ may indeed 
decrease very rapidly with an increase of mo for large values of the later. There are two 
reasons behind the above behavior: i) the choice of vanishing Higgs scalars in the NUSM 
at the unification scale Mq and ii) a large slepton mass (LSM) effect in the RGE of Eq.4 
via the tau- Yukawa term. Henceforth, we will refer the combined effect of i) and ii) as the 
LSM effect. In this analysis we first compute fnjj^ analytically without ignoring the terms 
involving hi, and hr- The calculation essentially keeps the terms involving the bottom and 
the tau- Yukawa couplings similar to what was used for the top- Yukawa term of Eq.4 in 
Refs. [6,35]. This results into, 

3 

m^^ = Cho^I + ml/^g{t) + -Sop. (10) 

Here g{t) is a function of tan/9 [35]. The term involving Sq vanishes in the NUSM. At the 
level of approximation where the bottom and tau Yukawa couplings are straightway ignored 
in Eq.4, would become zero in the NUSM. The result of our computation of Ch]j that 
leads to a non-vanishing value is given below. 

CHn=CHAU)+CHANU) (11) 

where, 

Ch^U) = 1-3(3/2 + 13) 

- h{K+Se,+SH^). (12) 

The quantities I2 and are functions of Yi = -^^^ for i = 2,3 (6, r). Small tan j3 solutions 
of YiS, computed at the electroweak scale [35] are shown in the appendix. I2 and /s defined 
below are computed numerically. 

h= fY2{t')dt\ h=fY^{t')dt' (13) 

JO JO 

Here t = ln{M^/Ml). One finds J2, /s << 1. For NUSM, the above reduces to 

Ch, = -2/3. (14) 
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Since 1^ is proportional to the square of r- Yukawa coupling where oc we find 
\Ch]j\ to be an increasing function of tan/5. For tan (3 = 10, we find — —0.007. The 
above clearly shows the large slepton mass effect because the dependence of Ch^ ''^'hd '^^ 
J3 arises from the term of Eq.4 that is associated with h"^. Thus with the NUSM parameters, 
the scalar mass term that appears as the first term in Eq.lO is essentially contributed by 
the third generation of slepton masses. Eq.lO and Eq.l4 clearly show that even with a small 
tan/3, larger values of rriQ may reduce m'jj^ appreciably and may turn the later negative. 
This demonstrates the LSM effect as mentioned before'^. We are not aware of any past 
reference that pointed out this large mass effect or provided with semi-analytic expressions. 
For larger tan /3, Cho ni£^y be appreciably large and negative that causes the LSM effect to 
become prominent even for a smaller value of niQ. 

We note that there is no LSM effect that may modify yU^ of Eq.7 at the one-loop level. 
This is simply because /i^ ~ —m^^^ and there is no possibility of having a LSM effect in 
the corresponding RGE of Eq.5. On the other hand, an improved rra^^ as obtained above 
in Eqs.lO and 14 changes the result of m\ (Eq.8 and 9) but the later does not receive any 
additional rn^ dependence other than what is already contributed from Eq.l4. Thus it follows 
that the coefficient of when Eq.8 is modified would be Ch]j- As a result, in the above 
analysis that is valid for small tan/? below 10 or 15, we find that the LSM effect causes m\ 
to have a HB/FP like behavior for its dependence on mo and mi. Hence, m^ may become 
significantly small for a large mo because of a cancellation between the terms. This of course 
may happen even for a small value of tan/3. A very large mo would cause m\ to become 
tachyonic, or would result into an absence of REWSB (Eq.3). For a large tan (3 on the other 
hand, the LSM effect is drastically enhanced. We comment here that in spite of showing the 
one-loop results we performed a complete numerical solution of the RGEs up to two loops 
in this analysis using SUSPECT [36]. 

We now point out that a simple non-universal Higgs scalar mass scenario as described 

^(i)We must mention that if we had chosen the values of Higgs scalar masses to be mo at the unification 
scale, the resulting coefficient of mg in Tn^j^ would be positive and the later would no longer be a decreasing 
function of jtiq- (ii) We further note that the fact that the third generation of squark masses are vanishing 
in NUSM also makes the LSM effect prominent. For example, had we considered non-vanishing values for 
mg and mjj {— tuq, ic. 6qj^ = S^i^ ~ 0) wc would find the terms in Eq.l2 that arc associated with I2 to 
become more prominent. This would have caused mj^^ to be further negative. However, it could also cause 
rriQ to turn tachyonic at the clectroweak scale. We have also checked this fact numerically. 
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in Ref. [10] may also provide a small rriA along with a funnel type of region that satisfies 
the WMAP data, for a smaller value of tan /5. We emphasize in particular the case where 
non-universality was analyzed with a unified Higgs scalar mass in Ref. [10]. However, unlike 
the NUSM these scenarios are very much constrained via REWSB so that considering larger 
values of masses for the first two-generation of scalars, an easier way to control FCNC and 
CP-violation effects is not possible. In these non-universal Higgs scalar models rriA becomes 
small for small values of tan f], only when tachyonic values of fn^Ho '^'hu (assumed equal) 
are considered at the unification scale and this severely reduces the available parameter space. 
Additionally, one may obtain an A-pole annihilation region or a funnel region for a small 
value of tan/3 in the so called sub-GUT CMSSM scenario [37]. 

We will now describe our results as obtained by using SUSPECT [36]. Fig.l shows the 
variation of /i and niA when mo is varied in the NUSM and mSUGRA. Fig. 1(a) shows the 
effect of varying mo in mSUGRA and in the NUSM for mi/2 = 500 GeV and 1 TeV. A 
small reduction of ((10 to 20%)) is seen in the NUSM when mo is increased up to 5 or 10 
TeV. The approximate one-loop result of Eq.7 and Eq.9 with 5 = — 1 as in NUSM however 
would indicate a flat /i over a variation of mo. The figure of course shows a moderately 
varying /i and we have checked that this variation has its origin in two-loop RGE effects. 
The two different mi/2 contours for mSUGRA however show a decreasing behavior of yU when 
mo is enhanced. This happens simply because of the HB/FP effect existing in mSUGRA. 
Fig. 1(b) shows the results for m^ which have two sets of contours for mi/2 = 500 GeV 
and 1 TeV corresponding to mSUGRA and the NUSM. In contrast to mSUGRA where mA 
rapidly increases with mo, the behavior in NUSM is opposite so that we may find a very light 
pseudoscalar Higgs boson via the LSM effect. The decrease of m^ is even more pronounced 
for a large value of tan/3 and this may be easily seen in Fig. 1(c) where a variation of tua 
vs mo is shown for a given mi/2(= 500 GeV) for three different values of tan/3. The curves 
ends in the larger mo sides for a few different reasons. For tan/3 ^ 10, the largest mo limit 
is caused by stop mass becoming very light or unphysical. On the other hand for tan /3 ^ 15 
the largest mg limit is given by absence of REWSB (ie m\ turning negative). 

The LSM effect may be explicitly seen in Fig. 2 where we plot squared scalar masses of 
Higgs scalars as well as their difference m\ ~ fnjj^ — m^^ with respect to a variation over 
the renormalization scale Q for mSUGRA and the NUSM for tan/3 = 10, mi/2 = 400 GeV 
and Aq = 0. Fig. 2(a) corresponds to mSUGRA for mo = 1 TeV whereas Fig.2(b) shows the 
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Figure 1: (a) Variation of fi with mo for tan j3 = 10, Aq = and mi = 500, lOOOGeV. The red 
dashed hnes represent the curves for the NUSM and the black dotted hues represent the same 
for mSUGRA. (b) Same as (a), but CP-odd Higgs (m^i) is plotted against mo- (c) shows 
the variation of with mo in the NUSM for different values of tan (3 when mi = SOOGeV. 

2 



case of the NUSM for tuq = 5 TeV. We have chosen different values of mg in the two figures 
because of the fact that a large mo is prohibited in mSUGRA via the REWSB constraint 
whereas a choice of a small mo in the NUSM would have a negligible LSM effect. Clearly, 
m'jj^ stays positive at the electroweak scale in mSUGRA (Fig.2(a)). On the other hand the 
same for the NUSM turns toward a negative value while running from Mq to the electroweak 
scale (Fig.2(b)) and this essentially shows the LSM effect in m\. Thus a typical large tan/3 
phenomenon that occurs in mSUGRA is obtained in the NUSM for a small tan (3. 
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Figure 2: (a) Variation of squared Higgs scalar masses and their difference (= m\ ~ m'jj^ — 
m?Hu) witli renormalization scale Q for tan/3 = 10, mi = 400 GeV, mo = 1000 GeV and 
Ao = for mSUGRA. (b) Same as (a), except for mo = 5000 GeV and for the NUSM . 



rrijj^ turns negative via LSM effect. 



3 Cold dark matter constraint and extended funnel re- 
gion 

In supergravity type of models Xi becomes the LSP for most of the parameter space [38,39] 
and we assume that the cold dark matter relic density is entirely due to Xi- Considering 
WMAP data [40] one finds a 3a limit as shown below. 



0.091 < VLcDMh'^ < 0.128 



(15) 



where QcDAih"^ is the DM relic density in units of the critical density. Here, h = 0.71 ±0.026 
is the Hubble constant in units of 100 Km s~^ Mpc~^. In the thermal description, the LSP 
was in thermal equilibrium with the annihilation products at a very high temperature of the 
early universe (T >> m^o). The annihilation products include fermion pairs, gauge boson 
pairs, Higgs boson pairs or gauge boson-Higgs boson combinations and they are produced via 
s, t and M-channel processes. As the temperature decreased the annihilation rate fell below 
the expansion rate of the universe and the LSP went away from the thermal equilibrium 
and freeze-out occured. The current value of ^xoh"^ is computed by solving the Boltzmann 
equation for n^o, the number density of the LSP in a Friedmann- Robert son- Walker uni- 
verse. The above computation essentially involves finding the thermally averaged quantity 
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< (JeffV >, where v is the relative velocity between two annihilating neutralinos and Ue// 
is the neutralino annihilation cross section that includes all the final states. In addition to 
annihilations one considers coannihilations [41-45] which are annihilations of the LSP with 
sparticles close in mass values with that of the LSP. The cross section sensitively depends 
on the nature of the composition of the LSP. In the MSSM, the LSP is a mixed state of bino 
(B), wino (W) and Higgsinos {H'l H2): 

X? = NnB + NuWs + iVigif? + iVi4^2 • (16) 

Here Nij are the elements of the matrix that diagonalizes the neutralino mass matrix. The 
gaugino fraction Fq of the lightest neutralino is defined by Fg = lA^'np + |A''i2p. A gaugino- 
like LSP may be defined to have Fg very close to 1(^ 0.9). On the other side, a Higgsino-like 
LSP would have Fg ^ 0.1. For values in between, the LSP could be identified as a gaugino- 
higgsino mixed state. The MSSM with gaugino masses universal at Mg has a few distinct 
regions in general that satisfy the WMAP constraint. In this section we point out the exis- 
tence of such regions in the context of the NUSM. i) The bulk annihilation region in mSUGRA 
is typically characterized by smaller scalar masses and smaller values of mi/2 where the LSP 
is bino-dominated. A bino dominated LSP couples favorably with right handed sleptons. 
Thus, the LSP pair annihilation in the bulk annihilation region occurs primarily via a t- 
channel sfermion in mSUGRA. There are two important constraints that disfavors the bulk 
region in mSUGRA. These are the constraints a) from the slepton mass lower limit from 
LEP2 [46] and b) from the lower limit of Higgs boson mass nih of 114.4 GeV [47]^. ii) The 
focus point [19] or the hyperbolic branch [20] region of mSUGRA is typically characterized 
by a small region that is close to the boundary of the lighter chargino mass lower bound. 
Because of a small \fi\ here the LSP has a significant amount of the Higgsino component 
or it can even be almost a pure Higgsino ^. Additionally the lighter chargino xt becomes 
lighter and coannihilations [43,44] with LSP reduce the relic density to an acceptable level. 
The HB/FP region is however absent in the NUSM. 

iii) Coannihilations of LSP may also occur with sleptons, typically staus (fi) [41] in mSUGRA. 
These regions are associated with small mi/2 and small mo zones near the boundary of the 
discarded zone where staus become the LSP. Stau coannihilation is also an effective way to 

^It is quite possible to have a bulk region with a non-zero Aq [48] . 
^This holds in the inversion region of the Hyperbolic branch [20]. 



13 



bring the neutralino relic density to an acceptable level in the NUSM. Coannihilations of 
LSP may also occur with stop (ti) in a general MSSM scenario [42] or even in mSUGRA [48]. 
However, in spite of having relatively lighter ti, the NUSM does not have such a region unless 
one reduces the mass of ti further via appropriately considering non-zero values for Aq. 
iv) The most important region satisfying WMAP data for our study is the Higgs-pole an- 
nihilation or funnel region [49,50]. The funnel region that satisfies the WMAP data is 
characterized by the direct-channel pole 2m^o ~ myi,m/^. This occurs in mSUGRA typi- 
cally for large tan/3 extending to larger mo and larger mi/2 regions. In the NUSM however, 
the funnel region occurs in all possible tan/9(^ 5). Indeed apart from the LSP-stau coan- 
nihilation appearing in a small region, Higgs-pole annihilation is the primary mechanism in 
the NUSM to satisfy the WMAP constraint throughout the parameter space. 

We now show the results of the computation of the neutralino relic density using the 
code micrOMEGAs [51]. In Fig.3(a) that corresponds to tan/? = 10 the region with red 
dots in the (mi — mg) plane satisfy the WMAP constraint (Eq.l5) for the neutralino relic 
density. For a small mg we find the stau coannihilation region marked with red dots. The 
gray region below the stau coannihilation region is discarded because of the appearance of 
charged LSPs. 

The upper gray region is discarded broadly via m^ turning negative except near the 
boundary where there can be interplay with other constraints as described below. Thus if 
we concentrate on the boundary of the discarded region, the smallest mi/2 zone (below 160 
GeV or so) is ruled out by the LEP2 lower limit of sparticle masses [46]. This is followed 
by obtaining tachyonic sfermion scalars (particularly stop scalars) when mi is increased 
further up to ^ 600 GeV. The same boundary zone for the larger mi/2 region is eliminated 
because of the appearance of the charge and color breaking (CCB) minima [52]. In the 
NUSM this happens via the CCB conditions that involve rn^Ho ' later becoming negative 
makes the CCB constraint stronger. In the region between the stau coannihilation area and 
the upper gray shaded discarded area one finds a long red region that satisfies the WMAP 
constraint for VL^oh? . As mentioned before, the reason for satisfying the WMAP data is the 
direct channel annihilation of two LSPs via neutral Higgs bosons. Fig. 3(b) shows a scanned 
output for tan/? = 10 when mo is varied as in Fig. 3(a). The LSP mass is plotted against 
(2mj^o — rriA) /2'm^o so as to show the extent of the A-annihilation or funnel effect. Clearly 
the WMAP satisfied points shown in red fall around the zero of the y-axis confirming that the 
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funnel region appears in the NUSM for small tan /3 as well. We note that the A-width can be 
quite large [Ta ~ 10-50 GeV) and 2m^o can be appreciably away from the exact resonance 
zone still giving a s-channel annihilation consistent with the WMAP data. The heavy scalar 
Higgs boson H also significantly contributes to the total annihilation cross-section. We now 
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Figure 3: (a) WMAP allowed region in the mi/2 — rriQ plane for tan/? = 10 and Aq = with 
yU > for NUSM. Lighter Higgs boson mass limits are represented by solid lines. Dot-dashed 
line refers to 6 ^ 57 limit. WMAP allowed regions are shown by red dots. The entire 
region is allowed via Bg f^'^fJ'" bound. The small figure within the inset is shown for an 
improved display of the stau coannihilation zone, (b) Scattered points in the plain of m^o 
vs (2m^o — m^)/2m^o shown after a scanning of mi/2 and mo for tan/? = 10. Almost all 
the WMAP satisfied points (in red) occur near the zero of the y-axis, thus suggesting the 
s-channel annihilation of the LSPs via A and H bosons. 



discuss the result of using a few constraints. First, we use the LEP2 limit [47] of 114.4 GeV 
for the SM Higgs boson. For the CP-even lighter SUSY Higgs boson mass rrih we use the 
same constraint as long as we are in the decoupling region [31] where the SUSY /i-boson is 
SM-like. In the NUSM this is true for almost all the parameter space except a very small 
region with small mi and for tan/3 ^ 15 which we will discuss later. Additionally, we 
note that there is an uncertainty of about 3 GeV in computing the mass of the light Higgs 
boson [53]. This theoretical uncertainty primarily originates from momentum-independent 
as well as momentum- dependent two-loop corrections, higher loop corrections from the top- 
stop sector etc. Hence, we have drawn a contour for rrih = 111 GeV in order to consider an 
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Figure 4: (a)WMAP allowed region in the mi/2 — niQ plane for tan/3 = 15 and = with 
/i > for NUSM. Lighter Higgs boson mass limits are represented by solid lines. Dot-dashed 
line refers to 6 ^ S7 limit. WMAP allowed regions are shown by red dots. The — ^ fi~^fi~ 
bound is represented by the long-dashed line and this discards a small strip of region below 
the discarded top gray region. The small figure within the inset is shown for an improved 
display of the stau coannihilation zone, (b) Variation of sin(/5 — a) vs mo for mi = 500 GeV, 
tan/3 = 15 and Aq = for NUSM. The figure shows that for large mo, actually near the 
region where niA is small, sin(/5 — a) is consistently small. This is of course much away from 
decoupling region of SUSY Higgs boson. The figure in inset shows the finer variation of 
sin(/3 — a) within a small range of mo. The gray shaded region refers to the similarly shaded 
discarded region of Fig. 4 (a). 



effective lower limit. ^ 

We have further drawn the fe — ^ S7 contour by considering a 3a limit [54] , 

2.77 X 10"^ < Br{b S7) < 4.33 x 10"^ 



(17) 



Clearly in Fig. 3 (a) this constraint would keep the mi/2 > 425 GeV region alive and this 
is indeed the region of super-large mo ^ 4 TeV that satisfies the WMAP data. In this 
context we should however note that the 6 — > s + 7 constraint may be of less importance 
if some additional theoretical assumptions are taken into consideration. The computation 
of the constraint in models like mSUGRA assumes a perfect alignment of the squark and 
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^The following values of top and bottom quarks are considered in the SUSPECT code used in our analysis: 
= 172.7 GeV and m6*^(mb) = 4.25 GeV. 
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quark mass-matrices. This essentially considers an unaltered set of mixing angle factors than 
the Cabibbo-Kobayashi-Maskawa (CKM) factors at the corresponding SM vertices. Even a 
small set of off-diagonal terms in squark mass matrices at the unification scale may cause 
a drastic change in the mixing pattern of the squark sector at the electroweak scale. This 
however does not cause any effective change in the sparticle mass spectra or in the flavor 
conserving process of neutralino annihilation or in generating events in a hadron collider in 
any significant way. A brief review of the model-dependent assumptions for the b— > S7 
analyses may be seen in Refs.( [55], [56]) and references therein. 

We now explain an interesting aspect of the NUSM Higgs boson as a consequence of the 
LSM effect. As we have seen before, rriA decreases with increasing mo in the NUSM. We 
find that with larger rriQ, ttia may become very light for the small mi zone so that we may 
find a spectrum where ~ ~ Ma and this indeed is a consequence of moving into 
the intense-coupling [32] region of Higgs bosons. This results in a change in the value of 
the coupling gzzh (oc sin(/9 — a)), where a is the mixing angle between two neutral Higgs 
bosons h and H. Typically in the decoupling region sin(/3 — a) stays very close to 1, but in 
the intense coupling region this can be considerably small like 0.5 or lesser. This effectively 
reduces the lower limit of mh to 93 GeV, a value close to mz- In the NUSM we obtain 
this effect for a small region of parameter space (small mi) if tan/5 ^ 15. Fig. 4(a) shows 
such a region (this may also satisfy the WMAP limits) for 200 GeV ^ mi ^ 300 GeV and 
this appears very close to the top gray shaded discarded zone. Fig.4(b) demonstrates the 
existence of a small sin(/9 — a) in the NUSM as discussed above. However with Aq = 0, 
we will see that such a very light or rrih region is almost discarded via the present limit 
of the Bg fi^fi^. The current experimental limit for the Br{Bs — > /i''^/U^) coming from 
CDF [57] puts a strong constraint on the MSSM parameter space. The experimental bound 
is given by (at 95 % C.L.) 

Br(fi, ^ /i+^-) < 5.8 X 10"^ (18) 

The estimate of Bs f^'^fJ'" [58] in the MSSM sensitively depends on the mass of A-boson 
(oc m^^) and on the value of tan (3 (oc tan^/3). Bs Z^"*"/^" constraint eliminates the thin 
region along the boundary of the REWSB in NUSM where m^ is very light. As mentioned 
before, for tan/5 ^ 15 a part of the above region for small mi zone has non-Standard Model 
like /i-boson. 

Apart from the above constraints we would like to remind that there is no non-universality 
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in the first two-generations in NUSM. This saves from the stringent FCNC violating hmits 
such as those coming from the Kl-Ks mass difference or from the /i — > 67 bound. Splitting 
of the first generation and the third generation of scalars may also cause violations of FCNC 
bounds, although to a lesser extent. Following Ref. [22] we see that for no violation of FCNC 
one would need (for equal gluino and average squark masses) \mq{l) — mg(3)| ^ m?/Miy. 
Here mq refers to the average squark mass. Considering the analysis performed in Ref. [22] 
for different gluino masses we conclude that the results of our analysis stay in the safe zone 
regarding the FCNC bounds in spite of the inter-generational splitting between the squarks. 
Fig. 5(a) shows the results for tan/3 = 40. Here the funnel region is extended up to mi/2 = 
1.7 TeV. A large value of tan (3 increases and hr and this would enhance the width Ta [59]. 
We point out that unlike Fig. 3(a), here the reach of mg decreases for a given mi/2- The mass 
of A-boson decreases with increase of mo much more rapidly because of larger and hj. 
arising out of larger values of tan (3. Here, the very light mA region that satisfies the WMAP 
data and that evades the LEP2 Higgs boson limit exists near the top gray boundary for 
250 GeV ^ nil ^ 350 GeV which is again ruled out via the Bs /U^yU" limit. The upper 
gray region is typically discarded here via the REWSB constraint of m\ that turns negative 
at the tree level. Finally, we have not imposed any limit from the muon g — 2 data that may 
or may not show a discrepancy from the Standard Model result. It is known that using the 
recent e^-e~ data leads to a 3.4(T level of discrepancy. On the other hand, using hadronic 
T-decay data in computing the leading order hadronic contribution to muon g — 2 washes 
away [60] any deviation from the SM result. 



NUSM 
tanp=40 
A„=(),|i>0 



b->s+7^ 







1000 
m,„ (GeV) 



a 

i_ 

a" -5- 



NUSM 

tanP=40,Ao=0,^>0 



400 600 
mo (GeV) 



(a) (b) 

Figure 5: Same as Fig. 3 except tan/3 = 40. 
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3.1 Sample parameter points satisfying WMAP data for early run 
of LHC 

We would focus on a few characteristic parameter points to discuss the nature of the NUSM 
spectra that satisfy the WMAP hmits. Clearly one is able to reach a considerably large mg 
satisfying dark matter relic density constraints and all other necessary constraints for mi < 
1 TeV. As an example with mi = 1 TeV, this limit is around 7 TeV for tan/5 = 10, 6 TeV 
for tan/5 = 15, and 1.6 TeV for tan/5 = 40. Thus, for larger tan/5 the reach of mo decreases 
considerably. On the other hand, it is not unusual to obtain an A-po\e annihilation region 
for a large tan/5 in popular models like mSUGRA. Hence, if we are interested in focusing 
on a region of the MSSM that is not available in mSUGRA type of models we would rather 
explore the smaller tan /5 domain of the NUSM. The other important zone of parameter space 
could be the region with smaller mi because this would be easily accessible in LHC in its 
early run. We pick up a point on Fig. 3 (a) (scenario-A) that just satisfies the minimum value 
of rrih = 111 GeV besides being consistent with the WMAP data for cold dark matter. The 
input parameters are: tan/5 = 10, Aq = 0, mi = 270 GeV, mo = 2050 GeV and sign{fi) = 1. 
The scenario-A of Table 1 thus has light stop and light sbottom quarks, light charginos 
and neutralinos and at the same time it would have a light Higgs spectra, all of which are 
promising for an early LHC detection. We remind ourselves that typically the NUSM is 
associated with a heavy first two-generation of scalars and heavy sleptons for all the three 
generations. We note that we have relaxed the 6 — > S7 constraint keeping in mind of the 
argument given after Eq.l7. We could of course respect the constraint, only at a price of 
having a little heavier spectra, still that would be very much accessible in LHC for the 
third generation of squarks, Higgs bosons, charginos etc. This however would not cause any 
essential change in the general pattern. On the other hand, with not so light m^ and with a 
small tan/5 the scenario-A satisfies the Bg — > fi'^fi~ limit. The scenario-B of Fig.4(a) refers 
to a special point (tan/5 = 15, Ao = 0, mi = 255 GeV, mo = 2000 GeV and sign{fi) = 1) 
for which the Higgs sector is not in the decoupling region thus reducing the limit of rrih to 
a value near the Z-boson mass. This parameter point obeys the Bg fi'^fi" limit but as 
with the scenario-A it has an inadequate Br{b 57). With smaller mi and smaller mo 
scenario-B has a lighter spectrum in general than scenario-A. Particularly, we find a much 
smaller Higgs boson spectra. The scenario-C (with tan/5 = 40, = 0, m 1 = 540 GeV, 
mo = 1250 GeV and sign{fi) = 1), that satisfy all the constraints including the b —>■ S'j 
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bound provides a relatively heavier spectra. This is however still within the LHC reach for 
the relevant part of the NUSM as mentioned above. We note that the large mo domain 
of the NUSM is typically associated with lighter Higgs spectrum and this region may in 
general be probed via the production of all MSSM Higgs bosons and subsequently their 
decays. The squarks corresponding to the first two generations and sleptons of all the 
generations are large when mo is large. As a result the production cross-section of these 
particles will be too low at the LHC. Thus, to see any SUSY signal for the NUSM one 
should primarily analyze productions and decays of gluino,stop and sbottom in addition 
to charginos and neutralinos. In a part of the parameter space where the Higgs spectra 
is light, all the Higgs states may be produced via pp /i, H, A,H^ + X either via loop- 
induced processes like {gg —>■ h, H, A) or through cascade decays via heavier charginos and 
neutralinos pp — > xt^ xh Xa ~^ Xi^ X2' Xi + H, A, H^. The latter decays are only allowed 
if enough phase space is available. 



4 Direct and Indirect detections of dark matter 
4.1 Direct detection rates 

We will now discuss the prospects of direct and indirect detections [61] of neutralino (LSP) as 
a candidate for dark matter in the NUSM . First, we will discuss the direct detection of LSP 
via measurements of nuclear recoil. Neutralinos interact via spin-independent (scalar) and 
spin-dependent interaction [38,62] with nucleons. The scalar cross-section may be expressed 
in terms of number of protons and neutrons, Z and [A — Z) respectively [38] as follows. 

(Tscalar = ^—[Zfp + (A - Z) Q\ (19) 
TT 

where is the reduced LSP mass. The quantities fp and /„ contain all the information of 
short distance physics and nuclear partonic strengths and these may be seen in Refs. [63]. 
We will now comment on the relative strengths of spin-independent and spin-dependent 
neutralino-nuclear cross-sections. While (Jscaiar depends on Z and A — Z quadratically, the 
spin-dependent interaction cross-section on the other hand is proportional to J{J + 1) where 
J is the total nuclear spin [38] . Typically the spin-independent neutralino- nucleon scattering 
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parameter 


A 


B 


C 


tan/3 


10.0 


15.0 


40.0 




270.0 


255.0 


540.0 


mo 


2050.0 


2000.0 


1250.0 


^0 











sgn{fi) 


1 


1 


1 




312.60 


291.53 


651.60 


m-g 


709.32 


674.43 


1278.26 




2103.25 


2047.68 


1658.64 




276.89 


248.27 


842.00 




493.53 


465.21 


1028.15 




390.66 


354.64 


958.59 




434.06 


403.08 


1019.74 




2050.19 


1999.70 


1296.33 




2037.46 


1972.60 


1119.32 




196.44 


183.99 


430.22 


"^xt 


347.34 


327.11 


668.46 




347.72 


326.97 


668.18 




318.11 


297.66 


655.13 


m~o 
X2 


197.69 


185.15 


430.30 


m~o 

Xl 


108.05 


101.72 


226.91 




259.48 


148.37 


403.03 




272.027 


169.44 


411.73 


rrih 


111.26 


111.25 


116.32 


^1 


0.105 


0.102 


0.13 




1.59 X 10"^ 


4.65 X 10-5 


2.73 X 10-^ 




4.02 X 10-9 


2.81 X 10-*^ 


5.29 X 10-« 




9.31 X 10-" 


1.59 X 10-1° 


6.99 X 10-1° 



Table 1: Data for point A, B, and C. Masses are in GeV 

cross-sections ( where cr^p^si — ^xn,si) are appreciably smaller than the corresponding spin- 
dependent cross-sections {ay^p^sD — o'xn,SD)- However considering the fact that aso oc J{J + 
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1) and a SI oc Z'^,{A — ZY, a scalar becomes considerably larger for moderately heavy elements 
{A > 30) [3,64] like Ge, Xe etc. But we should keep in mind that there exists some cases 
where asn may become considerably larger than asi even for A > 30. 

The cross-section ascaiar mainly involves the computation of x ~ Q and x ~ 9 scattering 
amplitudes. The scalar cross-section at tree level is composed of t— channel Higgs boson 
exchange and s— channel squark exchange contributions. On the other hand, the spin- 
dependent cross-section depends on t— channel Z exchange and s— channel squark exchange 
diagrams. In this analysis we compute the spin-independent cross-section agcaiarix ~ p) ^^i 
two values of tan /3 (10 and 40) for = and ;U > using DARKSUSY [65]. Fig. (6) shows 
o^scaiarix ~ p) ^s m^o for tan/? = 10 and 40 when mi/2 and mo are varied as in Fig. 3 and 
Fig. 5 respectively. The WMAP satisfied points are shown in small (maroon) circles. We have 
shown the limits from CDMS (Ge) 2005 [66], XENON-10 [67] and from future experiments 
like SuperCDMS (Snolab) [68] and XENONIT [69]. The neutralinos with higher mass (up to 
400 GeV) will be partially probed in XENON-IT while the light LSP region for tan p = AO 
is already ruled out by XENON-10 data. 
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Figure 6: Spin-independent scattering LSP-nucleon cross-sections vs LSP mass for the NUSM 
for tan /3 = 10 and 40. The blue dotted region correspond to scanned points of the parameter 
space of Figs. 3 and 5. The maroon filled circles show the WMAP allowed points. Various 
limit plots are shown for different experiments. 
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4.2 Indirect detection via photon signal 

The fact that the A-resonance annihilation is the primary mechanism to satisfy the WMAP 
hmits in the NUSM suggests that there will be enhanced signals for indirect detection via 
7-rays, positrons and anti-protons in the NUSM . On the other hand, detection via neutrino 
signal [61] would not be interesting in this case where the LSP is almost a bino annihilating 
via A-resonance in the s-channel [59]. Among the above indirect detection possibilities we 
will limit ourselves to estimating only the detection prospect of gamma rays that originates 
from the galactic center [70-73]. In general for neutralino annihilation at the galactic center 
one may have the following possibilities: i) monochromatic 7-rays and ii) continuum 7-rays. 
Monochromatic 7-rays come out from processes like XX ~* 11 [74] and XX ~^ ^1 [75]. These 
signals although small because of the processes being loop suppressed are clean with definite 
energies = and = — m|/4m^. Continuum 7-rays on the other hand arise from 
neutralinos annihilating into a variety of Standard Model particles. The hadronization and 
production of neutral pions would follow. Typically pion decay, in particular 7r° 77 would 
produce a huge number of photons with varying energies. 

The differential continuum 7-ray flux that arrives from angular direction ip with respect 
to the galactic center is given by [61,71,72], 



Here cxj is a LSP pair annihilation cross section into a final channel i. We will consider 
7-rays emerging from the galactic center, hence if) = v is the pair's relative velocity and 
< av > refers to the velocity averaged value of av. is the differential 7-ray yield for 
the channel i. p^i^) is the cold dark matter density at a distance r from the galactic center, 
where = + -Rg — 2sRqcosiP. Here, s is the line of sight coordinate, Rq is the Solar 
distance to the galactic center. Clearly, Pj;(r) that depends on astrophysical modelling is 
important to determine the photon flux in Eq. (20). One may indeed preferably isolate the 
right hand side of Eq. (20) into a part depending on particle physics and a part depending 
on astrophysics. For the later, one defines a dimensionless quantity J{ip) such that. 




(20) 



0.3 GeV/cm'^y June of sight 




dspi{r{s,^)). 



(21) 
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The above results in, 

— ^(E^, tb) = 0.94 X 10~^^cm"2s"^GeV ^sr"^ V — ^— — - J(th). 

dE^^ ^ rfE^ V10-29cm3s-V V ) 

(22) 

For a detector that has an angular acceptance Af2 and lowest energy threshold of Eth the 
total gamma ray flux from the galactic center is given by, 

2 



(23) 

Here J(Af2) = /^^^ J{ilj)dQ. The upper limit of the integral in Eq.(23) is flxed by the 
fact that the neutralinos move with galactic velocity, therefore the annihilations may be 
considered to have occured at rest. We will now comment on the galactic halo density 
proflles used in this analysis. Various N-body simulations suggest that one may obtain a 
general proflle behavior arbitrary to the extent of a few parameters and this is given by [76], 

pW^P. I^ + Wl"^..,, . (24) 
(r/iio)"'|l + (r/a)'>] — 

Here po is a normalisation factor which is taken as the local {i.e. solar region) halo density 
(~ 0.3 GeV/cm^). We will analyze with three popularly used proflles, namely the isothermal 
cored [77], Navarro, Frenk and White (NEW) proflle [78] and Moore proflle [79] as given in 
the Table (2). The table also mentions the corresponding value of J for AQ = 10""^ and 10~^ 
sr. Computation of the photon flux for a different halo proflle may easily be performed by an 
appropriate scaling with the corresponding J. Clearly, more cuspy proflles would produce 
higher photon-flux. One can further include the effects of baryons on the dark matter halo 
proflles. Baryons may undergo radiative processes leading to a fall towards the central region 
of a galaxy in formation. This changes the density proflles of matter towards the center which 
in turn leads to an increased concentration of dark matter. Adiabatic compression [80] has 
been used to study the baryonic effects. Inclusion of the adiabatic compression effects cause 
the proflles to become signiflcantly cuspier, often increasing J by a factor of 100 or so [72]. 
We have not included these halo proflle models in our computation, but the photon flux 
would increase by a similar factor as mentioned above. 
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Halo Model 


a (kpc) 


Ro (kpc) 


a 


/9 


7 


J(10-3) 


J(10-5) 


Isothermal cored 


3.5 


8.5 


2 


2 





30.35 


30.4 


NFW 


20.0 


8.0 


1 


3 


1 


1.21 X 10^ 


1.26 X 10^ 


Moore 


28.0 


8.0 


1.5 


3 


1.5 


1.05 X 10^ 


9.75 X 10^ 



Table 2: A few dark matter halo density profiles and associated parameters. 
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Figure 7: Continuous 7-ray flux in cm~^s~^ above a threshold energy of 1 GeV for a cone 
of 1 X 10^^ sr centered around the galactic center vs rriQ. Lines are shown for three different 
halo distributions i) spherically symmetric isothermal cored profile (SPH) [77], ii) Navarro, 
Frenk and White (NFW) profile [78] and iii) Moore profile [79]. 
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Fig. (7) shows the result of continuum photon- flux vs mo for tan/5 = 10 and 40 correspond- 
ing to two different values of mi/2 (= 400 and 800 GeV) and three different halo profiles as 
mentioned above. The photon flux (in cm~^s~^) in the NUSM is computed using DARK- 
SUSY [65] with > 1 GeV, for a solid angle aperture of AQ = 10~^ sr. As mo increases 
the mass of psuedoscalar Higgs boson decreases and the LSP pair annihilation through s- 
channel resonance causes a peak corresponding to a halo profile. Broadly, such a peak covers 
the region of mo for a given mi/2 where WMAP data for the neutralino relic density is 
satisfied. With an increase in tan (] the width of rriA increases and the peak associated 
with a given halo profile broadens. We see that in spite of having a broad range of halo 
profile characteristics, the photon flux in the region of resonance annihilation in the NUSM 
where the WMAP data is satisfied may be probed in the upcoming GLAST [81,82] exper- 
iment at least for the cuspier profiles. GLAST would be able to probe photon-fiux as low 
as 10^^° photons/cm^/s [82]. We note that this conclusion remains valid in spite of the fact 
that GLAST will use an aperture of AQ = 10"^ sr so that an appropriate scaling of the 
photon-fiux in Fig. (7) needs to be done from Table (2) and Eq.(23). We further note that 
as mentioned before, the use of the adiabatic compression mechanism would modify a given 
halo profile significantly and this may increase the photon flux by a few orders of magnitude. 
Fig. (8) shows the plots of photon-fiux (in cm~^s~^ with > 1 GeV) vs the mass of the 
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Figure 8: Scatter plot of photon flux in cm~^s~^ (with E^ > 1 GeV) vs LSP mass for 
tan/3 = 10 and 40 for a NFW halo profile in the NUSM . Here mi/2 and mo are scanned in 
the ranges shown in Figs. (3,5). Only WMAP relic density satisfied points are shown. 



LSP for tan/? = 10 and 40 for NFW halo profile in the NUSM . Here, mi/2 and mg are 
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varied such that mi/2 < 2 TeV and mo < 20 TeV for tan/5 = 10 and mi/2 < 2 TeV and 
mo < 10 TeV for tan/3 = 40. Only WMAP allowed parameter points are shown. GLAST 
would be able to probe up to 400 to 450 GeV. 

5 Conclusion 

In this analysis we have worked with a non-universal scalar mass scenario in a supergravity 
framework. We started with purely phenomenological motivations namely, i) to manage the 
FCNC and CP-violation type of constraints by decoupling, ii) to obtain WMAP satisfied 
values for neutralino relic density for a broad region of parameter space without depending 
on any delicate mixing of gauginos and Higgsinos, iii) to have radiative electroweak sym- 
metry breaking and iv) to keep naturalness within control. Keeping the above in mind and 
considering a unified gaugino mass scenario, we used a common scalar mass parameter mo 
at the gauge coupling unification scale for the first two-generation of scalars as well as the 
third generation of sleptons. The item (i) mentioned above would require mo to be large, 
and the item (iv) would prefer light third generation of squarks and light Higgs scalars. For 
simplicity we used vanishing third generation of squarks and Higgs scalar masses at the uni- 
fication scale. In such a scenario with a possibly multi-TeV mo, we first found a large mass 
effect or in particular large slepton mass effect in the ROE of m^^ [ie. for the down type 
of Higgs scalar) that turns the later negative at the electroweak scale almost irrespective of 
a value of tan/3. Extending the semi-analytic solution of m\ by considering the hh and the 
hr terms (this is required here even for a small tan/3) relevant for the large slepton mass 
effect we found a hyperbolic branch/focus point like effect in m^ (^ ""^/f^ — ^"h^j) for small 
tan/3. This causes ttla to be almost independent of mo for a large domain of the later. But, 
with a very large mo this causes rriA to become very light or this may even turn m\ negative 
giving rise to no radiative electroweak symmetry breaking. We further found that because 
of such large slepton mass effect in the RGE, the Higgs sector may reach an intense coupling 
region with all the Higgs bosons becoming very light and this may evade the LEP2 limit 
of rrih for a limited region of parameter space. However, constraint from Br{Bs — > /U"'"/i~) 
becomes stringent in this region because of lighter m^. In general, we find relatively lighter 
rriA OT rriH and this fact leads to a large A-pole annihilation region or funnel region of dark 
matter even for a small tan/3. This is in contrast to minimal supergravity type of scenar- 
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ios where funnel region may occur only for large values of tan [3. The nature of the LSP 
is bino-dominated, thus there is no need of any delicate mixing of binos and Higgsinos in 
order to satisfy the neutralino relic density constraint. We have also computed the direct 
detection rates of LSP-nucleon scattering. The upcoming detectors like XENON-IT would 
be able to probe almost the entire region of parameter space. We have further estimated the 
indirect detection prospect via computing continuous photon fluxes. The ongoing GLAST 
experiment will be succesfully able to probe the parameter space even for a less cuspy halo 
profile. We have also briefly discussed the detection prospect of sparticles in the LHC The 
Higgs bosons and the third generation of squarks are light in this scenario. In addition to 
charginos and neutralinos the above may be easily probed in the early runs of LHC. 



6 Appendix 

The coefficients appearing in Eqs. 7 and 8 are given by, 



^ tan^/? , ^ 1 , ^ 2m 1^ tan^/5 _ 

^^= (tan^/3-l) ^' ^^ = - (tan^/3-l) ^^-^^"" = " (tan^ /3 - 1) ^ ^^S) 

tan^/j + 1 tan2/5 + l tan^/^ + l 

Here the functions k,g,e and / may be seen in Ref. [35]. Electroweak scale result of 
Yi = hf/^ATi)^ with i = t,b, and r are shown below. 

m)YM EM) ..dFsW^l^sW (27) 
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The quantities Ei{t) are defined as follows. 

16 3 13 

Er{t) = (l + /33(t))^(l + /32(t))^(l + /5i(t))^ 

E3(t) = (l + /3i(t))^(l + /32(t))* (28) 

Here, F{t) = Ei{t')dt', A = tti(0)6i/47r and (6i, 62, ^3) = (33/5, 1,-3). 
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